MMCTrem.  calculations  of  the  response  of  aspheru.„™. 

PAVLOAD  TO  CONING  NOTIONS)  RRNV  BALLISTIC  RESEARCH  LAB 
ABERDEEN  PROVING  GROUND  ND  G  R  COOPER  HAR  88 
UNCLASSIFIED  BRL-MR-JG61  F/G  19/18  NL 


BRL-MR-366  I 


Wc  FIIE 


MEMORANDUM  REPORT  BRL-MR-3661 


BRL 


1938  -  Serving  the  Army  for  Fifty  Years  -  1988 


NUMERICAL  CALCULATIONS  OF  THE  RESPONSE 
OF  A  SPINNING  SPHEROIDAL  PAYLOAD  TO 
CONING  MOTION 


GENE  R.  COOPER 


DTIC 

^  ELECTEB 
JUN  2  2  1988 1 

*  *h  1 


MARCH  1988 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


U.S.  ARMY  LABORATORY  COMMAND 

BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN  PROVING  GROUND,  MARYLAND 


88 


Kt 


15  6 


DESTRUCTION  NOTICE 


Destroy  this  report  when  it  is  no  longer  needed.  DO  NOT  return  it  to  the 
originator. 


Additional  copies  of  this  report  may  be  obtained  from  the  National  Technical 
Information  Service,  U.S.  Department  of  Ccrmerce,  Springfield,  VA  22161. 


The  findings  of  this  report  are  not  to  be  construed  as  an  official  Department 
of  the  Army  position,  unless  so  designated  by  other  authorized  documents. 


The  use  of  trade  names  or  manufacturers'  names  in  this  report  does  not  con¬ 
stitute  indorsement  of  any  ccmnercial  product. 


la.  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2*.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION  /DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUM8£R(S) 

BRL-MR-3661 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


Form  Approved 
OMB  NO.  0704-01  $9 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

U . S .  Army 

Ballistic  Research  Laboratory 


6c  ADDRESS  (City.  State,  and  ZIP  Code) 

Aberdeen  Proving  Ground,  MD 

21005-5066 

8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION  u<s<  Anny 

Ballistic  Research  Laboratory 

8b.  OFFICE  SYMBOL 
(If  applicable) 

SLCBR-DD-T 

8c.  ADDRESS  (City,  State,  and  ZIP  Code) 

Aberdeen  Proving  Ground,  MD 

21005-5066 

3.  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release; 
distribution  is  unlimited. 


5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


7*.  NAME  OF  MONITORING  ORGANIZATION 


7b.  ADDRESS  (City.  State,  *nd  ZIP  Cod') 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO 

62618A 


62618AH80 


1 1  TITLE  (Include  Security  Classification) 

NUMERICAL  CALCULATIONS  OF  THE  RESPONSE  OF  A  SPINNING  SPHEROIDAL  PAYLOAD  TO  CONING  MOTION 


12.  PERSONAL  AUTHOR(S) 

Cooper,  Gene  R. 


13a.  TYPE  OF  REPORT 

Memorandum  Report 


16.  SUPPLEMENTARY  NOTATION 


13b  TIME  COVERED 
FROM _ TO 


14.  DATE  OF  REPORT  (Year,  Month.  Day) 

1988  March 


FIELD 

GROUP 

01 

02 

COSATi  CODES 


SUB-GROUP 


18  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

'Singular  Integral  Rotating  Liquid 

Numerical  Calculations,  Eigenvalue 

Spheroidal  Container  _ 


19.  ABSTRACT  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 


A  method  for  numerically  evaluating  a  particular  type  of  singular  integral  coupled  with 
a  complex  eigenvalue  problem  is  presented  here.  This  presentation  provides  the  information 
needed  to  carry  out  the  numerical  calculations  found  in  the  cited  Murphy  report.  This  tech¬ 
nique  has  proven  useful  for  calculating  the  moment  exerted  by  a  low-viscosity  spinning  liquid 
in  a  coning  spheroidal  container.  .  , 


20  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 
□  UNCLASSIFIED/UNUVITED  E  SAME  AS  RPT 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Gene  R.  Cooper 


121  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED  _ 


22b  TELEPHONE  (Include  Area  Code)  22c.  OFFICE  SYMBOL 

(301) -278-2057  SLCBR-LF-A 


t»4»lJUX4LUM^:LTamiJ!CT 


Table  of  Contents 


f 


Page 


I.  INTRODUCTION 


II.  NUMERICAL  INTEGRATION 


III.  ROOT  FINDER 


IV.  DISCUSSION 


REFERENCES 


LIST  OF  SYMBOLS 


DISTRIBUTION  LIST 


Copy 

CSP[CT£0 


Accession  For 


:  \  ^  i  i  bu 1 1  cn/  _ 

/.  -J-.*  ■  i-'il-  ilil-y  C o 3 e> 3 

i  and/or 
•  .-.I  ,  j;.*clal 


I.  INTRODUCTION 


Murphy1  has  studied  the  moment  exerted  by  a  low- viscosity  spinning  liquid  in  a  coning 
spheroidal  container  and  presented  some  computations.  The  purpose  of  this  report  is  to 
elaborate  on  some  of  the  numerical  methods  used  by  Murphy.  He  indicates  the  need  of 
further  clarification  of  two  numerical  problems.  The  first  problem  involves  the  numerical 
evaluation  of  integrals  that  have  a  singularity  in  their  integrands.  The  second  problem 
deals  with  calculating  complex  zeroes  (i.e., complex  eigenvalues)  of  a  complex  function. 


II.  NUMERICAL  INTEGRATION 


The  calculations  used  in  Reference  1  basically  center  around  numerically  evaluating 
integrals.  These  integrals  are  found  in  Equations  (6.22)-(6.23),  (8.6)  and  (8.9)  of  his  report 
and  are  repeated  here  for  convenience: 


(Ukm(skm)),  +  Re-l/2akm  =  0 


(6.22) 
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Equation  (6.22)  represents  the  normal  velocity  boundary  condition  that  must  be  sat¬ 
isfied  at  the  wall  of  the  container  for  non-forced  motion.  ( £4m),-  is  the  (A:,m)th  mode  of 
the  normal  inviscid  velocity  for  k  =  2,4,6  and  m  —  0, 1  .  The  variable  a*m  is  obtained 
by  expanding  the  normal  wall  boundary  layer  correction  in  the  least  squares  formulation 
given  by  the  terms  Ckkm  and  bkkm  ■  Expression  (8.6)  gives  the  non-dimensional  side  moment 
due  to  the  inviscid  pressure  plus  the  boundary  layer  correction  to  this  pressure  resulting 
from  the  velocity  boundary  layer  at  the  wall.  The  boundary  layer  also  introduces  a  fur¬ 
ther  correction  to  the  liquid  side  moment  found  from  the  viscous  diffusion  terms.  This 
correction  is  given  by  Equation  (8.9).  The  definitions  as  well  as  the  interpretation  of  the 
remaining  variables  in  the  above  equations  can  be  found  in  Murphy’s  report  and  the  List 
of  Symbols  at  the  end  of  this  report.  The  corresponding  integrations  are  carried  out  over 
the  range  of  [—1,1]  .  By  inspection  it  is  revealed  that  each  integrand  is  an  even  function 
of  u>  .  Hence  each  integration  is  equal  to  twice  that  integration  over  the  interval  [0, 1]  .  For 
constant-amplitude  motion,  e*m  =  0,  Murphy  points  out  that  the  integrands  in  Equations 
(6.23)  and  (8.6)  have  a  singularity  at  u>  =  ±u>c  ,  where 


(2.1) 


uc  =  /  |  m  -  r  |  [4  +  /a(m  -  rfe]  1/2 
m  =  0, 1 

After  some  straightforward  algebra  it  can  be  shown  that  each  integral  in  Equations 
(6.23)  and  (8.6)  can  be  written  in  the  form: 

j  =  r  n»,°„0,)d»  +  /■  F(u...<i-'-i)ai,,Kift)ju, 

./o  "s/^c  —  U  iut  \/ U)  —  <*)c 

where  F(w,  qi,  ft)  is  an  analytic  function  on  [0, 1].  Equation  (2.2)  shows  explicitly  the  type 
of  singularity  found  in  the  corresponding  integrands.  Identify  the  first  integral  in  Equation 
(2.2)  by  II  and  second  by  12.  Using  the  substitution  t2  =  uc  —  ui  in  II  and  t2  =  u  —  uc  in 
12  allows  the  integral  I  to  be  written  as: 

rZZ  f\/ 1— u»e 

I  =  2i( — i)C1— lmD  /  G(t,  ■  ■  -)dt  +  2  /  G(t,---)dt  (2.3) 

*/o  */o 

The  integrals  shown  in  Equation  (2.3)  are  no  longer  singular.  These  integrals  along  with  the 
integral  given  in  Equation  (8.9)  can  now  be  integrated  by  standard  numerical  quadratures. 
The  particular  numerical  quadrature  that  was  used  for  Reference  1  is  an  iterative  Simpson’s 
rule  with  adaptive  gridding.  This  is  a  variation  of  a  technique  given  by  Sampine  and  Allen2 
where  the  real  integrand  is  replaced  by  a  complex  integrand. 


III.  ROOT  FINDER 

The  results  given  in  TABLE  3  of  Reference  1  can  now  be  found  by  using  a  numerical 
root-finding  technique  on  Equation  (6.22).  For  completeness  TABLE  3  is  reproduced  on 
page  4. 

The  procedures  described  above  allow  the  values  of  this  analytic  complex  function  to 
be  calculated.  The  eigenvalues  (Murphy’s  Sfcm)  were  readily  found  by  applying  a  Muller3 
algorithm  to  Equation  (6.22).  Its  main  advantage  is  that  it  requires  functional  values  only, 
not  derivatives.  Experience  has  shown  that  with  reasonable  initial  guesses  for  smk,  the 
Muller  algorithm  obtains  the  desired  roots  with  only  a  few  iterations. 

IV.  DISCUSSION 


The  procedures  presented  in  this  report  show  one  method  of  eliminating  the  singular 
point  in  the  required  integrals.  These  integrals  as  well  as  the  eigenvalues  can  then  be  eval¬ 
uated  using  standard  methods  found  in  text  books  on  numerical  analysis.  The  numerical 
methods  presented  here  were  found  to  work  successfully  when  applied  to  the  Murphy  1 
integration  and  the  eigenvalue  problems.  All  of  the  remaining  calculated  results  found  in 

2  Samptne  and  Allen,  1973, Numerical  Computing: A n  Introduction  ,W,B.  Saundert  Company. 

3Mulier,  H,,  1956:  A  Method  for  Solving  Algebraic  Equation*  Using  an  Automatic  Computer.  Math.  Comput.  10,208-215. 


Table  1:  Comparison  of  Eigenvalues  with  Greenspan/Kudlick  (TABLE  3.  of  ref.  1) 
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